NPS-OR-92-013 


AD-A256  073 


NAVAL 


POSTGRADUATE  SCHOOL 

Monterey,  California 


An  Application  of  Growth  Curve  Analysis  to  the 
Ammunition  Stockpile  Deterioration  Model 


SO  YOUNG  SOHN 


August  1992 


Approved  for  public  release;  distribution  is  unlimited 


Prepared  for: 

Naval  Weapons  Support  Center 
Crane,  IN  47522 


92-26359 

in  mu  mu  mu  mu  mu  mu  m:  hi 


NAVAL  POSTGRADUATE  SCHOOL 
MONTEREY,  CALIFORNIA 


Rear  Admiral  R.  W.  West,  Jr.  Harrison  Shull 

Superintendent  Provost 

This  report  was  prepared  for  and  funded  by  the  Naval  Weapons 
Support  Center,  Crane,  Indiana.  Reproduction  of  all  or  part  of 
this  report  is  authorized. 

This  report  was  prepared  by: 


SO  YOUNG  SOHN 
Assistant  Professor 
Operations  Research  Department 


PETER  PURDUE 

Chairman  and  Professor 

Operations  Research  Department 


PA 
Dean  of 


search 


UNCLASSIFIED 

security  class >.cat'0\  this  page 


REPORT  DOCUMENTATION  PAGE 


la  REPORT  SECUR-y  C.aSSiF.Cat  ON 


id  RESTRICTIVE  MARKINGS 


2a  SECURITY  C.ASS  F  CATION  AuThOR. ty 


3  distribution  availability  of  report 


2b  DEC.ASSif iCAHON  •  DOWNGRAD  NG  SCHEDULE 


4  PERFORMING  ORGANIZATION  REPORT  NUMBER(S) 


8 


proved  for  public  release;  distribution  unl- 


5  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 


S-OR-92-013 


6a  NAME  OF  PERFORMING  ORGANIZATION  6b  OFFICE  SYMBOL 

(If  applicable ) 

Naval  Postgraduate  School  OR 


6c.  ADDRESS  (Ofy,  State,  and  ZIP  Code) 

Monterey,  CA  93943-5000 


7a  NAME  OF  MONITORING  ORGANIZATION 


Naval  Postgraduate  School 


7b  ADDRESS  (C/ty,  State,  and  ZIP  Code) 

Monterey,  CA  93943-5000 


8a  NAME  Oc  E  jNDNG  SPONSORING 
organization  Naval 
pons  Support  Center 


8c.  ADDRESS  (City,  State,  and  ZIP  Code) 

Zrane ,  IN  47522 


8b  OFFICE  SYM80L 
(If  applicable) 


9  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 

RRJKA 


io  source  of  funding  numbers 


PROGRAM 
ELEMENT  NO 


PROJECT 

NO 


11  Title  (Include  Security  Classification ) 

Sn  Application  of  Growth  Curve  Analysis  to  the  Ammunition  Stockpile  Deterioration  Model 


12  PERSONAL  AuThOR(S) 


>  •  W  f  >f  1 


13a  type  Oc  REPORT 


1 3b  time  covered 

:,0MJul  92  T0  Sep  92 


14  DATE  OF  REPORT  (Year,  Month  Day)  115  PAGE  COUNT 


COSA'I  CODES 


GROUP 


SJB-GROUP 


18  SUBJECT  TERMS  ( Continue  on  reverie  if  necessary  and  identify  by  block  number) 

Random  Effect  Logistic  Regression  Model,  Deterioration  Rate, 
A  Two-Stage  Estimation 


19  ABSTRACT  ( Continue  on  reve'se  if  necessary  and  identify  by  block  number) 

Ammuniation  deterioration  during  storage  has  considerable  economic  consequences.  A  reliable 
prediction  model  for  the  ammunition  deterioration  rate  is  necessary  for  long-term  procure¬ 
ment  and  maintenance  planning.  A  random  effect  growth  curve  analysis  is  employed  to  formulat 
a  predition  model  for  amnunition  deterioration  rates  in  terms  of  concurrent  characteristics 
such  as  depot  condition  and  vendor  information.  The  resultant  prediction  model  can  be  used 
to  determine  the  appropriate  time  for  reorder  or  renovation  of  anmunition  before  ths  quality 
reaches  unacceptable  levels.  A  two-stage  analysis  is  used  to  estimate  parameters  involved 
in  the  prediction  model.  Necessary  estimation  methods  are  discussed.  An  example  is  given 
to  illustrate  the  imp lenient at ion  procedure  of  the  prediction  model  suggested  in  this  paper. 


20  D'ST R'BUTiON  AVAILABILITY  OF  ABS'RACt 

jj]  UNCLASSIFIED  UNLIMITED  □  SAME  AS  RPT  □  DTlC  USERS 


22a  NAVE  0s  RE $PONS;Blt  iND  VDLAL 


21  ABSTRACT  SECUR'TY  CLASSIFICATION 


l»J  i  ’•#  F-tCrt  I 


22b  TELEPHONE  (Include  Area  Code)  22 c  OFFICE  SvMBOL 


tc'sm'/vi  ii  '.t.  •c'si  m 


DD  FORM  1473,  84  va  = 


83  APR  ed  !  on  "nay  be  used  until  exhausted 
A  !  other  ecitior-s  are  obsolete 


SECUR'TY  classification  of  this  page 

O  j  S  Government  Printing  0**»ce  1fl( — 


Unclassified 


An  Application  of  Growth  Curve  Analysis 
to  the  Ammunition  Stockpile 
Deterioration  Model 


So  Young  Sohn 

Department  of  Operations  Research 
Naval  Postgraduate  School 
Monterey,  CA  93943 


Abstract 

Ammunition  deterioration  during  storage  has  considerable  eco¬ 
nomic  consequences.  A  reliable  prediction  model  for  the  ammunition 
deterioration  rate  is  necessary  for  long-term  procurement  and  main¬ 
tenance  planning.  A  random  effect  growth  curve  analysis  is  employed 
to  formulate  a  prediction  model  for  ammunition  deterioration  rates  in 
terms  of  concurrent  characteristics  such  as  depot  condition  and  vendor 
information.  The  resultant  prediction  model  can  be  used  to  determine 
the  appropriate  time  for  reorder  or  renovation  of  ammunition  before 
the  quality  reaches  unacceptable  levels.  A  two-stage  analysis  is  used 
to  estimate  parameters  involved  in  the  prediction  model.  Necessary 
estimation  methods  are  discussed.  An  example  is  given  to  illustrate 
the  implementation  procedure  of  ;he  prediction  model  suggested  in 
this  paper. 

Key  Words:  Random  Effect  Logistic  Regression  Model,  Deterioration 
Rate,  A  Two-Stage  Estimation 
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1.  INTRODUCTION 

Each  year,  the  U.S.  military  services  purchase  millions  of  dollars  worth 
of  major  caliber  ammunition  from  several  vendors  to  prepare  for  possible 
war.  Ammunition  lots  purchased  are  stored  in  depots  until  they  are  used.  In 
view  of  the  fact  that  most  ammunition  is  produced  long  before  its  ultimate 
consumption,  it  is  important  that  the  material  be  adequately  stored  and 
maintained  to  remain  usable. 

However,  according  to  the  GAO  report  (The  Comptroller  General  (1982)), 
the  military  services  are  experiencing  significant  problems  in  maintenance 
of  depot  and  renovation  of  ammunition.  These  problems  cause  ammunition 
to  deteriorate  faster  than  necessary  maintenance  can  be  performed  and  a 
serious  backlog  of  ammunition  in  need  of  maintenance.  In  an  attempt  to 
prevent  such  undesirable  phenomena,  in  this  paper,  a  quantitative  model 
that  can  relate  a  depot  maintenance  plan  to  the  quality  control  of  ammu¬ 
nition  is  considered. 

The  objective  of  this  paper  is  to  outline  a  prediction  model  for  ammuni¬ 
tion  deterioration  in  terms  of  depot  characteristics  and  vendor  information. 
The  goal  is  to  provide  estimation  methods  by  which  one  can  determine  the 
appropriate  time  for  reorder  or  renovation  of  ammunition  before  the  qual¬ 
ity  reaches  unacceptable  levels.  The  suggested  model  can  also  be  used  for 
vendor  control  and  examination  of  depot  maintenance  policy. 

In  order  to  formulate  such  a  prediction  model  for  ammunition  dete¬ 
rioration  rates,  a  random  effect  growth  curve  analysis  is  employed.  The 
analysis  of  a  random  effect  growth  curve  model  has  been  extensively  used 
in  biomedical  research  (Rao  (1965),  Stiratelii  et  al.  (1984),  and  Vonesh 
and  Carter  (1987)).  An  application  of  growth  curve  analysis  to  physical 
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engineering  is  given  in  Sohn  and  Mazumdar  (1991).  In  general,  the  goal  of 
growth  curve  analysis  is  to  estimate  random  growth  rates  of  a  time  varying 
principal  attribute  in  terms  of  concurrent  individual  characteristics  which 
do  not  change  over  the  experimental  period.  The  analysis  of  growth  curves 
is  typically  carried  out  by  using  a  two-stage  random  effect  model.  The 
first  stage  consists  of  the  within-individual  regression  model  in  which  the 
within-individual  growth  rates  for  the  serial  observations  are  estimated  for 
all  individuals.  In  the  second  stage,  estimated  growth  rates  are  related  to 
a  set  of  covariates  representing  concurrent  individual  characteristics. 

In  this  paper,  a  growing  pattern  of  defective  proportion  of  an  attribute 
of  ammunition  lot  over  time  (deterioration  rate)  is  viewed  in  analogy  to 
a  growth  rate  of  a  principal  attribute  of  an  individual  experimental  unit 
while  vendor  information  and  depot  characteristics  associated  with  each 
lot  can  be  considered  as  concurrent  individual  characteristics.  In  order  to 


predict  deterioration  rates  of  a  certain  type  of  ammunition  in  terms  of  depot 
characteristics  and  vendor  information,  a  random  effect  logistic  regression 
model  (Stiratelli  et  al.  (1984))  is  employed  assuming  that  the  pattern  of 
deterioration  follows  a  logistic  function.  For  estimation  of  parameters  in 
a  random  effect  logistic  regression  model,  a  two-stage  estimation  based  on 
the  method  used  in  Korn  and  Whittemore  (1979)  is  employed.  An  example 
is  given  to  illustrate  the  implementation  procedure  of  the  random  effect 


logistic  model  suggested. 
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2.  RELATED  STUDIES 

The  most  important  deterioration  processes  occurring  to  ammunition 
during  storage  are  corrosion  of  the  metallic  components  and  chemical  and 
physical  changes  within  the  internal  structure  of  the  explosives,  which  may 
result  in  reduced  stability  and  performance  (Eckman  (1974),  Eriksen  and 
Stromspe  (1978),  and  Wilson  and  Bird  (1986)).  The  factors  which  have  the 
greatest  influence  on  these  processes  are  temperature  and  humidity  which 
are  directly  related  to  the  location  and  condition  of  the  depot.  Even  under 
similar  conditions,  different  types  of  ammunition  may  degrade  with  differ¬ 
ent  deterioration  rates.  In  some  instances,  different  vendors  provide  the 
same  types  of  ammunition  that,  nonetheless,  have  varying  deterioration 
rates.  Most  of  the  studies  in  the  literature  do  not  reveal  the  quantita¬ 
tive  relationship  between  deterioration  rates  and  depot  condition,  types  of 
ammunition  and  vendor  information. 

For  instance,  Eriksen  and  Strpmsoe  (1978)  conducted  an  extensive  ex¬ 
periment  covering  a  9  year  period  to  examine  deterioration  of  six  types  of 
artillery  ammunition  (105  mm  HE,  90  mm  HE,  57  mm  HE,  40  mm  HE, 
hand  grenades  Mk2,  and  antitank  mines  M6A2)  at  three  aboveground  and 
three  underground  depot  sites.  Various  attributes  were  observed:  vacu¬ 
um  stability  and  water  content  of  propellant,  impact  sensitivity  of  primers, 
burning  rate  of  propellant,  muzzle  velocity,  and  the  number  of  misfires  and 
duds.  Inspection  frequencies  varied  from  two  to  four  times  covering  com¬ 
binations  of  starting  year  (1968)  and  the  subsequent  1st,  5th,  7th  and  9th 
year  after  storage,  depending  upon  the  attributes  inspected.  Based  on  these 
repeated  measurements,  the  contingency  table  and  the  analysis  of  variance 
were  used  to  test  if  the  levels  of  ammunition  deterioration  were  significantly 
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different  in  terms  of  the  conditions  of  depot  and  the  types  of  ammunition. 

The  results  of  their  study  are  as  follows:  Overall,  aboveground  sites 
surrounded  by  inland  climate  performed  better  than  underground  depots 
that  were  old  and  technically  in  poor  shape;  the  visual  inspection  revealed  a 
considerable  difference  in  corrosion  attack  between  the  ammunition  types; 
water  content  both  in  the  primer  and  the  propellant  charges  varied  signif¬ 
icantly  depending  on  the  depot  in  which  lots  of  ammunition  were  stored; 
most  of  the  chemical  analyses  were  of  little  interest  and  should  be  postponed 
to  the  end  of  trials.  Their  approach,  however,  did  not  provide  a  compre¬ 
hensive  model  which  related  the  level  of  deterioration  of  certain  type  of 
ammunition  to  depot  condit:on.  If  the  level  of  deterioration  over  time  can 
be  expressed  as  a  function  of  depot  conditions  or /and  vendor  information, 
reorder  time  of  a  certain  type  of  ammunition  which  reaches  certain  quality 
can  be  predicted  in  terms  of  given  information  at  the  depot.  An  appli¬ 
cation  of  growth  curve  analysis  which  facilitates  such  a  prognostication  is 
illustrated  in  the  following  section. 
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3.  A  GROWTH  CURVE  ANALYSIS  FOR  AMMUNITION 
DETERIORATION  RATE 

Consider  a  certain  type  of  ammunition  (say,  105mm  HE).  In  this  section, 
the  assumptions  used  in  the  growth  curve  analysis  to  predict  the  deteriora¬ 
tion  rate  of  105mm  HE  ammunition  are  summarized.  Similar  assumptions 
can  be  made  for  other  types  of  ammunition. 

Basic  unit  of  purchase  is  a  lot  and  as  a  result  of  acceptance  sampling, 
the  qualities  of  incoming  ammunition  lots  would  be  homogeneous  regardless 
of  different  sources  of  manufacturers.  While  several  hundreds  of  such  lots 
are  stored  in  various  locations,  the  quality  of  ammunition  deteriorates  with 
the  passage  of  time.  Deterioration  rates  may  vary  depending  upon  both  the 
condition  of  depot  in  which  ammunition  lots  are  stored  and  manufacturers 
who  provide  the  ammunition  lots. 

In  order  to  inspect  the  pattern  of  deterioration  of  various  lots  over  the 
time,  it  is  assumed  that  a  set  of  initially  selected  lots  in  various  depot- 
s  is  marked  and  these  lots  are  repeatedly  inspected  on  a  sampling  basis 
without  rectification  over  a  certain  time  period.  Lot  size  is  assumed  to 
be  sufficiently  large  that  it  would  not  be  depleted  at  the  end  of  repeated 
sampling  inspections.  Modes  of  inspection  can  cover  visual  inspection,  lab¬ 
oratory  examination  and  functional  trials  (Eriksen  and  Stromsoe  (1978)). 
The  proportion  of  samples  which  do  not  satisfy  a  predetermined  level  of  an 
attribute  or  the  combinations  of  attributes  can  be  used  as  a  proxy  for  the 
accumulated  level  of  deterioration  at  the  time  of  inspection.  For  instance, 
among  many  other  attributes  of  functional  trials,  the  relative  number  of 
misfires  and  duds  to  the  sample  size  observed  at  each  inspection  cycle  can 
be  considered  as  a  measure  for  the  level  of  deterioration  of  each  lot  at  the 
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time  of  inspection. 

For  each  lot,  the  proportion  of  defective  rounds  is  assumed  to  increase 
as  time  goes  on  following  a  logistic  curve.  The  logistic  model  is  well  es¬ 
tablished  as  a  basis  for  analyzing  such  a  phenomenon  when  the  response 
variable  follows  binomial  distribution.  The  logistic  model  formed  to  de¬ 
scribe  the  deteriorating  pattern  for  each  lot  is  called  a  within-individual 
model.  Deterioration  rate  of  a  lot  of  ammunition,  which  represents  a  dete¬ 
riorating  pattern  over  time  of  a  lot  often  differs  from  that  of  another  lot. 
The  model  in  which  varying  deterioration  rates  between  lots  are  compared 
in  relation  to  depot  characteristics  and  vendor  information  associated  with 
each  lot  is  called  the  between-individual  regression  model.  Depot  condition- 
s  and  vendor  information  associated  with  each  lot  of  certain  ammunition 
are  assumed  to  remain  constant  during  the  experimental  period  while  they 
would  vary  over  different  lots.  Necessary  notations  for  a  two-stage  random 
effect  logistic  model  for  ammunition  deterioration  rate  is  as  follows: 

NOTATIONS 

i:  index  for  individual  lot,  i  =  1, ..,  N 

N:  total  number  of  ammunition  lots  employed  in  the  experiment 

j:  index  for  inspection,  j  =  1, ..,  n, 

nt:  the  number  of  inspections  to  be  done  on  a  lot  i  over  experimental 
period 

jth  inspection  time  of  a  lot  i 

Y^:  sample  size  of  lot  i  at  the  jth  inspection 

x,jk:  a  dichotomous  variable  indicating  if  each  item  k  taken  from  the 
sample  at  the  jth  inspection  of  lot  i  is  defective  or  not 

y 

Xij :  the  number  of  defectives  found  in  a  sample  size  of  1^;  i  x«jA.- 
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pv  =  Xij/  Yij-.  proportion  defective  at  the  jth  inspection  of  a  lot  i 
exp(0o)/(l+exp(/3<)))-  initial  proportion  defective  of  ammunition  lot 
0i :  the  deterioration  rate  of  lot  i  (log  odds  ratio  over  the  unit  time 
interval  for  lot  i) 

(z«i,  Zim):  covariates  indicating  concurrent  characteristics  of  lot  i 
(e.g.,  dummy  variables  indicating  different  combinations  of  depot  condition 
and  vendor  information  of  ammunition  lot  i) 

7i, 7m:  regression  coefficients  for  (jzn, zim ) 

m:  total  number  of  covariates  in  the  between-individual  model 

7:  a  pxl  vector  of  72 ’s,  (7i,--,7m)' 

Z:  an  Nxm  matrix  of  covariates  (2,1, .. ,2,m ) 

€,:  random  error  of  logfl,  ~  N(0,  <r2). 

p estimated  expected  proportion  defective  of  lot  i  at  time  t  given 

(■^fl  1  -m  Zin,') 

f, :  estimated  expected  time  that  quality  of  lot  i  reaches  a  predetermined 
level  of  proportion  defective  p  given  (2,-], zim) 

f3o,  (3,:  maximum  likelihood  estimators  of  (30  and  /3, 
log/3:  an  iVxl  matrix  of  logj3x' s,  ( logi3\,..,logdu )' 

6,:  estimation  error  associated  with  log(3t,  6,  ~  N(0,  7 f) 

7^,  r,2:  variance  of  (30  and  that  of  0, 

7^,  f2:  maximum  likelihood  estimators  of  and  r2 

W:  an  NxN  diagonal  matrix  with  diagonal  elements,  1  / (o'2  +  r* /0])'s 

7:  maximum  likelihood  estimator  for  7 

Mathematical  formulations  regarding  the  random  effect  logistic  model 
are  as  follows: 

In  the  within-individual  model,  proportion  defective  of  lot  i  at  the  jth 


inspection,  p*,  is  assumed  to  follow  a  logistic  function  of  time  t$\ 

For  i=l, N,  and  j  =  1, n, 

Pij  =  exp{J3o+  +  exp(j30  +  (3^))  (1) 

or 

l°9(Pij/(\  -  Pij))  =  Po  +  Pity,  (2) 

In  the  betweeu-individual  model,  deterioration  rate  of  ammunition  lot 
2,  Pi ,  is  represented  as  a  function  of  covariates  (e.g.,  depot  condition,  and 
vendor  information)  as  well  as  a  random  error  (e,): 

Pi  =  f  {hi 2  a  +  ••  +  7 rn2im  +  fi)  (3) 

where  follows  independent  N(0,  a2).  The  proportion  defective  of  ammu¬ 
nition  lot  would  not  decrease  as  time  goes  on.  The  following  function  is 
used  to  ensure  positive  deterioration  rates,  i.e.,  positive  pi’s: 

P,  =  exp(y,Zi\  +  ..  +  hmZ,m  +  Ci)  (4) 

or 

log(pi)  =  2,-1  +  ..  +  7 mzim  +  e,  (5) 

After  replacing  Pt  in  (2)  with  that  in  (4),  one  could  estimate  regression 
parameters  ( p0 ,  71,..,7m,  and  cr2),  by  maximizing  the  marginal  likelihood 
of  the  data  based  on  the  combined  model.  However,  exact  solutions  for  the 
normal  equations  based  on  the  combined  model  are  analytically  and  com¬ 
putationally  infeasible  (Stiratelli  et  al.  (1984)).  Indicating  these  problems, 
Stiratelli  et  al.  (1984)  disc  -sed  alternative  ways  of  estimating  parameters 
(EM  algorithm  and  a  two  stage  estimation).  When  Pi  in  (2)  is  replaced  with 
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that  in  (4),  random  error  e/s  appear  as  coefficients  of  the  regression  model. 
In  EM  algorithm  they  are  treated  as  missing  observations  and  are  esti¬ 
mated  by  cycling  back  and  forth  between  the  E-step  and  the  M-step  until 
convergence  is  reached.  When  a  two-stage  estimation  method  is  used,  the 
estimation  procedure  for  the  within-individual  model  is  separated  from  that 
for  the  between-individual  model.  For  instance,  in  the  within-individual  re¬ 
gression  model,  deterioration  rates  of  a  set  of  individual  ammunition  lots 
are  estimated  for  all  experimental  lots.  In  the  between-individual  regression 
model,  the  estimated  deterioration  rates  are  used  as  values  of  a  dependent 
variable  and  they  are  regressed  on  depot  characteristics  and  vendor  infor¬ 
mation  associated  with  each  individual  lot. 

In  general,  a  two-stage  estimation  method  is  computationally  less  inten¬ 
sive  and  analytically  easier  than  the  EM  algorithm  (Stiratelli  et  al.  (1984)). 
In  the  following  section,  a  two-stage  estimation  is  described  in  order  to  for¬ 
mulate  a  prediction  model  for  ammunition  deterioration. 
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4.  A  TWO-STAGE  ESTIMATION 


First,  in  order  to  estimate  deterioration  rates  in  the  within- individual 
model  (1),  the  following  likelihood  function  of  x^k  conditional  on  0o,  0u--,0N 
is  formulated: 

i=nnn ^ -  n  n  ( ?  Wo - ?«)*- x ■  m 

*=lj=lk=l  i=l;=l  \a9/ 


where  p,y  is  as  in  equation  (1). 

Maximum  likelihood  (ML)  estimates,  /30,  0i,..,0n  axe  obtained  by  solv¬ 
ing  a  set  of  normal  equations  of  log  likelihood  function  of  (6)  with  respect 
to  0O,  Estimated  parameters  (/3,’s)  are  in  turn  assumed  to  be  ap¬ 

proximately  nc.mal  with  means  equal  to  the  true  individual  parameters 
(fii  s)  and  variances  (rf’s).  Diagonal  elements  of  the  inverse  of  the  nega¬ 
tive  information  matrix  evaluated  at  (flo,  0i,..,0n)  provide  the  estimated 
variances  of  0O,  0\ which  are  (t^,  t f , t£). 

Secondly,  resultant  0,  replaces  unobservable  0i  in  the  between-individual 
model.  This  replacement,  however,  brings  the  estimation  error  <5,  to  the 
equation  (5): 


log(/3,)  =  7l  Zfi  +  ..  +  ImZim  +  e.-  +  (7) 

where  <5,  is  statistically  independent  of  e,.  It  is  assumed  that  <5,  follows  in¬ 
dependent  Normal  distribution  with  mean  0  and  variance  if  /0]  (var(<5,)  = 
x?a(log0i  |  00  ~  var ((/?,— /?,)//?,)  =  if /02i).  The  var(tf,-)  is  approximated  by 
Ti/0i  and  log0,  is  marginally  distributed  as  independent  N(7iZ,]-h  ..,+7 mzim, 
a1  -I-  f i/0i).  Standard  normal  theory  can  be  used  to  estimate  the  between- 
individual  model  parameters  ji's  and  cr2  which  maximize  the  following  log 
likelihood  function  at  the  estimated  if  : 
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N  N  ^  m 

LL  =  -^2 1°9(^  +  *? )  -  J2(log(3t  -  £  nuzuf/io-2  +  ij/ff).  (8) 

*=i  »=i  i=i 

The  resultant  MLE  of  7  is 

7  =(Z'WZ)-\Z'W{logP))  (9) 

and  the  estimated  variance  of  7  is  {Z'WZ)~l  where  an  NxN  diagonal 
matrix,  W,  consists  of  diagonal  elements,  l/(d2  +  7f//?2)’s.  In  order  to 
be  able  to  estimate  7,  the  number  of  lots  ( N )  employed  in  the  experiment 
must  be  lager  than  the  number  of  covariates  (m)  to  be  used  in  the  between- 
individual  model.  In  addition,  as  N  increases,  accuracy  related  to  the 
inference  concerning  7  would  increase. 

When  these  tj's  replace  Vs  in  (7),  logfii  can  be  predicted  in  terms  of 

Zil  ,  ..,  -Zfm- 

log(4i)  =  71-2.1  +  ••  +  ImZim-  (10) 

Use  of  (10)  to  quality  control  of  ammunition,  vendor  control  and  depot 
maintenance  plan  is  illustrated  in  the  next  section  along  with  a  numerical 
example. 
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5.  IMPLIMENTATION 

Unfortunately,  actual  data  related  to  type,  quantity,  and  depot  condi¬ 
tion  of  ammunition  were  omitted  from  the  GAO  report  (The  Comptroller 
General  (1982))  to  keep  it  unclassified.  In  the  Norwegian  study  cited  in  the 
earlier  part  of  this  paper,  repeated  measurements  on  the  deterioration  of 
ammunition  was  not  exposed  for  all  the  experimental  lots.  In  fact,  only  the 
summary  statistics  for  deterioration  of  several  groups  of  ammunition  were 
available,  respectively.  In  order  to  illustrate  implementation  procedures  for 
the  methods  suggested  in  this  paper,  a  numerical  example  is  given  based  on 
the  simulated  data  generated  partly  based  on  the  parameters  used  in  the 
Norwegian  study  (Eriksen  and  Stromsoe  (1974))  and  the  guideline  of  the 
U.S.  Army  Ammunition  Surveillance  Procedure  (Supply  Bulletin  SB  742-1 
(1988)). 

Table  1  contains  information  regarding  the  series  of  the  number  of  defec¬ 
tive  items  ( Yq)  found  in  20  ammunition  lots  ( N ).  These  lots  are  classified 
into  four  groups  depending  upon  the  vendor  who  provided  ammunition  and 
depot  in  which  ammunition  is  stored:  group  1  consists  of  5  lots  (i  =  1, ..,  5) 
which  are  provided  by  vendor  A  and  stored  in  depot  1;  group  2  consists  of 
5  lots  (i  =  5, ..,  10)  which  are  provided  by  vendor  B  and  stored  in  depot  1; 
group  3  consists  of  5  lots  ( i  =  11,. .,15)  which  are  provided  by  vendor  A 
and  stored  in  depot  2;  and  group  4  consists  of  5  lots  (i  =  16,  ..,20)  which 
are  provided  by  vendor  B  and  stored  in  depot  2.  All  the  lots  are  inspected 
annually  based  on  sample  size  of  20  (X;)-  Total  number  of  inspections  of 
each  lot  (nt)  varies  over  4  different  groups:  n,  =  7  for  i  =  1, ..,  5;  n,  =  8  for 
i  =  6, ..,  10;  n,  =  9  for  i  =  11, ..,  15;  and  n*  =  9  for  i  =  16, ..,  20. 

Our  goals  are  first  to  find  the  relationship  between  deterioration  rates 
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and  cone  rrrent  characteristics  of  ammunition  lot  (depot  condition  and  ven¬ 
dor  information);  secondly  to  determine  if  any  group  of  ammunition  causes 
significantly  different  deterioration  rates  from  the  others;  and  finally  to  pre¬ 
dict  the  time  for  renovation  of  ammunition  in  terms  of  characteristics  of 
depot  and  vendor. 

In  order  to  achieve  the  first  goal,  a  two-stage  estimation  is  employed. 
First  of  all,  the  ML  estimates  for  parameters  in  the  within-individual  model 
0o,  and  f02,  t2,..,?^,)  are  obtained  based  on  (6).  For  this  step, 

PROC  LOGISTIC  of  a  statistical  package  SAS  (SAS  (1989))  is  used  by 
applying  the  values  of  ATy  ,  X),  Uj  in  Table  2  to  model  (2).  Resultant  /3's 
and  t2's  are  summarized  in  Table  2  while  /30  turns  out  to  be  -4.9034  which 
indicates  that  the  estimated  initial  proportion  defective  of  ammunition, 
(exp0 0) / (1  +  exp0 0))),  is  about  0.74%  regardless  of  sources  of  vendor. 
Sample  patterns  of  the  actual  deterioration  (Actual: X^j  Y \j  verses  t0)  are 
given  in  Figure  1  along  with  those  of  the  estimated  deterioration  (Fitted: 
(exp0o  +  4<<o)/(  1  +  exp0o  +  /Mr,))  verses  t0. 

The  between-individual  model  is  formed  to  relate  deterioration  rates  to 
the  ammunition  characteristics,  logfli  is  used  as  dependent  variable  and  the 
four  dummy  variables  (ztl,..,z;4)  representing  4  combinations  of  depot  and 
vendor  characteristics  are  used  as  covariates  without  an  intercept:  z,i  =  1 
for  ammunition  lots  which  belong  to  group  1,  za  =  0  otherwise;  za  =  1  for 
ammunition  lots  in  group  2,  =  0  otherwise;  z&  —  1  for  ammunition  lots 

in  group  3,  z^  =  0  otherwise;  2*  =  1  for  ammunition  lots  in  group  4,  z*  =  0 
otherwise.  Given  such  (z,i,  ..,2*)  as  well  as  (/?',  and  f2' s),  an  optimization 
package  GAMS  (GAMS,  1988)  is  applied  to  (8)  in  order  to  obtain  71, ..,  74 
and  <x2.  Estimates  Eire  summarized  in  Table  2  along  with  standard  errors  of 
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7i, 74  obtained  from  the  square  root  of  the  diagonal  elements  of  ( Z '  WZ)  1 

/ 15  0  0  0  \ 

0  15  0  0 

where  a  20x4  covariate  matrix  Z  —  ;  Is  =(1,1, 1,1,1)’  is 

0  0  15  0 

\  0  0  0  15/ 

a  5x1  unit  vector;  and  W  is  a  20x20  diagonal  matrix  consisting  of  diagonal 
elements  l/(<r2  +  f, •//??).  Values  of  if  and  $f's  are  given  in  Table  1  while 
a2=0.012. 

Estimated  7  (/  =  1,  ..,4)  is  the  weighted  average  of  logfi'fi  which  belong 
to  the  group  l.  A  large  7  indicates  fast  deterioration.  For  instance,  the 
estimated  average  log  deterioration  rate  of  ammunition  lots  in  group  2  (- 
0.8336)  is  lager  than  that  of  ammunition  lots  in  group  1  (-0.9022).  To  test 
if  any  group  characteristic  causes  significantly  faster  deterioration  rate  than 
another,  a  (1  -  a)  100%  confidence  interval  for  Z'7  can  be  used: 


where  l  is  a  4x1  vector  which  contrasts  a  pair  of  two  group  characteristics; 
t  is  a  student  t  distribution  with  degrees  of  freedom  N  —  p. 

In  order  to  test  if  the  difference  in  two  expected  log  deterioration  rates 
between  group  1  and  group  2  is  significant,  one  can  use  V  =  (1,  —1,0,0)  to 
obtain  a  95%  confidence  interval  for  71  —  72: 

(1,  -1,0, 0)(— 0.9022,  -0.8336,  -0.7184,  -0.6567)' 

± t( 0.025, 20  -  4)  v/(l,-l,0,0)(Z'lVZ)-1(l, -1,0,0)' 

/  0.00708  0  0  0  \ 

0  0.00473  0  0 

where  (Z'WZ)~l  = 

0  0  0.00278  0 

\  0  0  0  0.00226/ 
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The  resultant  interval  (—0.0686  ±2.12(0.10870))  contains  0  which  im¬ 
plies  that  at  q  =  5%  vendor  effect  is  not  significant  on  the  deterioration  of 
ammunition  stored  in  depot  1.  Similar  comparisons  can  be  made  by  mod¬ 
ifying  I's.  Complete  analyses  of  95%  confidence  intervals  for  other  pairs 
are  summarized  in  Table  3.  Overall,  significant  differences  are  not  found 
except  for  the  difference  between  group  1  and  group  4. 

The  resultant  estimates  from  (10)  can  also  be  applied  to  the  analyses  of 
ammunition  lots  which  are  not  used  in  the  experiment.  In  order  to  predict 
the  expected  quality  of  a  new  ammunition  lot  i  at  time  t  in  terms  of  depot 
and  vendor  characteristics  associated  with  lot  i  the  following  p )  is  used: 

p\  =  exp0o  +  i3,t)/(  1  +  exp0o  +  0{t))  (12) 

where  f3{  =  exp{jiZf  i+, ,  ±7  mz,’m).  For  instance,  the  predicted  proportion 
defective  of  ammunition  lot  i  provided  by  vendor  A  after  5-year  storage  in 
depot  1  is 

p )  ~  exp{~ 4.9034+ea:p(— 0.9022)x5)/(l+ezp(— 4.9034+exp(— 0.9022)x5)) 

=0.0534.  (13) 

In  addition,  one  can  estimate  the  expected  time  t }  when  the  quality  of 
ammunition  lot  {  reaches  a  predetermined  level  p: 

=  [log(p/(  1  -  p))  -  y3o]/exp(7iZ,-j+, ,  +7 mziTn).  (14) 

Expected  time  fe  when  the  quality  of  ammunition  lot  ?  provided  by 
vendor  A  and  stored  in  depot  1  reaches  p=  0.5  is  estimated  as 

t)  =  [log( 0.5/(l  -  0.5))  +  4.9034] /exp( -0.9022)  =  12.1.  (15) 
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For  the  four  groups  considered  in  this  example,  p)  when  5  and  f*  for 
p  —  0.5  are  given  in  Table  4. 
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6.  SUMMARY  AND  DISCUSSION 

A  random  effect  logistic  regression  model  was  used  to  formulate  a  predic¬ 
tion  model  for  ammunition  deterioration  rate  in  terms  of  depot  and  vendor 
information.  It  was  assumed  that  once  lots  are  selected  for  initial  inspec¬ 
tion,  they  will  be  inspected  repeatedly  over  a  certain  time  period  without 
any  rectification.  The  proportion  defective  was  considered  a  measure  of 
ammunition  deterioration.  Regardless  of  different  sources  of  vendors,  the 
initial  quality  of  ammunition  was  assumed  to  be  homogeneous  as  was  its 
eventual  deterioration  with  the  passage  of  time.' "The  log  odds  ratio  due  to 
unit  time  passage  was  defined  as  a  deterioration  rate  which  determines  the 
shape  of  the  logistic  curve.  Deterioration  rates  of  different  lots  of  ammuni¬ 
tion  often  vary  and  their  variation  was  assumed  to  be  explained  partly  by 
a  function  of  depot  characteristics  and  vendor  information  associated  with 
each  lot. 

A  random  effect  logistic  regression  model  was  used  to  predict  the  dete¬ 
rioration  rate  of  ammunition  in  terms  of  concurrent  characteristics  of  the 
ammunition  lot.  A  two-stage  analysis  was  employed  to  estimate  deterio¬ 
ration  rates  and  relate  them  to  the  characteristics  involved  in  ammunition 
lot.  Using  the  approach  of  this  paper,  one  can  determine  the  appropri¬ 
ate  time  for  reorder  or  renovation  of  ammunition  before  the  quality  reaches 
unacceptable  levels  when  depot  and  vendor  information  associated  with  am¬ 
munition  lot  is  available.  By  comparing  the  expected  deterioration  rates 
of  ammunition  stored  in  different  locations  one  can  improve  depot  mainte¬ 
nance  policy.  The  similar  approach  can  be  used  for  vendor  control. 

In  the  two-stage  model  used  in  this  paper,  we  considered  a  certain  type 
of  ammunition  and  its  deterioration  was  compared  in  terms  of  vendor  and 
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depot  characteristics.  In  order  to  compare  deterioration  rates  of  different 
types  of  ammunition,  the  following  approach  can  be  used:  First,  use  a  group 
of  relatively  homogeneous  types  of  ammunition  in  the  within-individual 
model,  and  secondly,  classify  them  not  only  in  terms  of  depot  condition 
and  vendor  information  but  also  by  types  of  ammunition. 

In  order  to  increase  the  accuracy  of  resultant  estimators  of  a  two-stage 
estimation,  sufficiently  large  number  of  inspection  which  captures  the  no¬ 
ticeable  deterioration  of  ammunition  would  be  necessary  over  many  exper¬ 
imental  lots.  However,  such  data  may  not  be  readily  available  in  reality. 
When  a  sufficient  number  of  repeated  measurements  of  actual  data  is  not 
available,  accelerated  testing  may  be  an  alternative  means  of  meeting  the 
time  and  budgetary  restrictions  of  the  experiment.  For  instance,  envi¬ 
ronmental  tests  under  simulated  climatic  extremes  are  an  integral  part  of 
normal  engineering  tests.  These  tests  include  accelerated  depot  conditions 
at  the  test  site  when  the  characteristics  of  the  test  item  necessitate  the 
gathering  of  such  data.  Nonetheless,  when  only  a  small  number  of  data 
is  available,  resampling  procedure  such  as  bootstrap  method  (Dalai  et  al 
(1990))  can  be  used  to  make  inferences  concerning  deterioration  rates  at 
the  cost  of  heavy  computation. 

The  measure  of  rates  of  ammunition  deterioration  considered  in  this 
paper  was  the  proportion  of  defective  items.  When  other  variables  such  as 
muzzle  velocity  or  water  content  are  used  as  a  measure  of  deterioration, 
modifications  of  the  within-individual  model  may  be  necessary  to  take  into 
account  the  different  nature  of  the  measure.  However,  the  basic  concept  of 
a  two-stage  analysis  still  applies  to  any  modified  models. 

A  restriction  used  in  the  between-individual  model  was  that  a  lot  of 
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ammunition  would  be  stored  in  the  same  depot  in  the  experimental  period. 
A  prediction  model  that  reflects  changes  in  depot  conditions  associated  with 
a  lot  of  ammunition  is  under  investigation  and  is  left  for  further  research. 


20 


REFERENCES 

Brooke,  A.;  Kendrick,  D.  and  Meeraus,  A.  (1988).  GAMS:  A  User’s 
Guide,  The  Scientific  Press,  Redwood  City,  CA. 

The  Comptroller  General  (1982).  "DOD  has  Serious  Problems  with 
Care  and  Maintenance  of  Conventional  Ammunition".  AD  A1 10764,  Re¬ 
port  to  the  Chairman,  Subcommittee  on  Defense,  Committee  on  Appropri¬ 
ations,  House  of  Representatives.  United  States  General  Accounting  Office, 
Washington  D.  C. 

Dalai,  S.  R.,  Fowlkes,  E.  B.  and  Hoadley,  B.  (1990).  "Risk  Analysis 
of  the  Space  Shuttle:  Pre-Challenger  Prediction  of  Failure".  Journal  of 
American  Statistical  Association,  84,  No.  408,  pp.  945-967. 

Eckman,  D.  R.  (1974).  "Economic  Comparison  of  Wood-Preservative 
Treated  and  Untreated  105mm  Ammunition  Boxes".  AD  A001532,  Army 
Armament  Command,  Rock  Island,  IL. 

Eriksen  S.  and  Stromspe  E.  (1978).  "An  Examination  of  Deterioration 
of  Ammunition  by  Storage,  Norwegian  Defense  Research  Establishment". 
AD  A055897,  Norwegian  Defence  Research  Establishment,  Norway. 

Korn,  E.  and  Whittemore,  A.  (1979).  "Methods  for  Analyzing  Panel 
Studies  of  Acute  Health  Effects".  Biometrics  35,  pp.  795-802. 

Rao,  C.  R.  (1965).  "The  Theory  of  Least  Squares  when  the  Parame¬ 
ters  are  Stochastic  and  its  Application  to  the  Analysis  of  Growth  Curves". 
Biometrika  52,  pp.  447-458. 

SAS  (1990).  SAS/STAT  User’s  Guide,  Vol.  2,  Version  6.  Fourth  Ed. 

Sohn,  S.  Y.  and  Mazumdar,  M.  (1991).  "A  Monte  Carlo  Study  of  Vari¬ 
able  Selection  and  Inferences  in  a  Two  Stage  Random  Coefficient  Linear 
Regression  Model  with  Unbalanced  Data",  to  appear  in  Communications 


21 


in  Statistics,  Theory  and  Methods  A20. 

Stiratelli,  R.,  Laird,  N.  and  Ware,  J.  H.,  (1984)  "Random-Effects  Models 
for  Serial  Observations  with  Binary  Response".  Biometrics  40,  pp.  961- 
971. 

Supply  Bulletin  SB  742-1,  Inspection  of  Supplies  and  Equipment:  US 
Army  Ammunition  Surveillance  Procedures,  Headquarters,  Department  of 
the  Army,  Washington  DC. 

Vonesh,  E.  F.  and  Carter,  R.  L.  (1987).  "Efficient  Inference  for  Random- 
coefficient  Growth  Curve  Models  with  Unbalanced  Data".  Biometrics  43. 
pp.  617-628. 

Wilson,  M.  A.  and  Bird,  R.  (1986).  "Defect  Investigation  of  Cartridge 
Illuminating  l-l/2inch".  AD  A173717,  Technical  Note,  Department  of  De¬ 
fense,  Defense  Science  and  Technology  Organization  Materials  Research 
Laboratories,  Melbourne,  Victoria,  Australia. 


22 


ACKNOWLEDGMENT 

I  am  grateful  to  Mr.  Larry  Massa  from  Naval  Weapon  Support  Center, 
Professors  Judith  Lind,  Keebum  Kang  and  Pat  Jacobs  from  the  Naval  Post¬ 
graduate  School  for  providing  me  with  a  great  deal  of  practical  information 
regarding  ammunition  deterioration  and  computational  methodology.  I  al¬ 
so  thank  Ms.  Hania  Laborn  for  her  editorial  assistance.  This  research  was 
partially  sponsored  by  U.S.  Naval  Weapon  Systems  Center,  Crane  IN.  The 
opinions  expressed  are  solely  those  of  the  author. 


23 


Proportion  Defective 


Actual 
O  Lot  1 
□  Lot  6 
A  Lot  11 
x  Lot  16 


Fitted 


Loti 
Lot  6 
Lot  11 
Lot  16 


Figure  1.  Deteriorating  Patterns  of  Ammunition  Quality 


Table  1:  Number  of  Defectives  X ,y  Observed  in  Sample  Size  of  20  at  the 
j'th  Inspection  of  Lot  i _ 


group  1  group  2 

lot  i 


tv 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

2 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

3 

1 

1 

1 

0 

1 

1 

1 

1 

1 

0 

4 

1 

1 

1 

1 

1 

1 

1 

1 

1 

0 

5 

1 

1 

1 

1 

1 

1 

1 

2 

1 

0 

6 

1 

2 

2 

1 

1 

2 

1 

3 

2 

1 

7 

2 

2 

2 

2 

1 

3 

2 

4 

3 

1 

8 

na 

na 

na 

na 

na 

3 

2 

5 

4 

1 

Pi 

0.3853 

0.4490 

0.4172 

0.3671 

0.3853 

0.4446 

0.3649 

0.4969 

0.4446 

0.2310 

% 

0.0071 

0.0055 

0.0063 

0.0077 

0.0071 

0.0033 

0.0045 

0.0028 

0.0033 

0.0088 

group  3 

lot  i 

group  4 

trj 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

2 

0 

0 

0 

0 

1 

1 

0 

0 

0 

0 

3 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

4 

1 

1 

1 

1 

2 

1 

1 

1 

1 

1 

5 

2 

1 

2 

1 

3 

2 

2 

2 

2 

r% 

£ 

6 

3 

2 

2 

2 

5 

4 

3 

3 

4 

2 

7 

4 

3 

3 

2 

7 

6 

4 

4 

5 

4 

8 

5 

4 

4 

2 

10 

8 

6 

6 

8 

5 

9 

7 

5 

6 

5 

13 

10 

8 

8 

10 

7 

10 

na 

na 

na 

na 

na 

13 

10 

10 

13 

9 

Pi 

0.4884 

0.4355 

0.4535 

0.4143 

0.6290 

0.5582 

0.5009 

0.5009 

0.5536 

0.4776 

% 

0.0019 

0.0022 

0.0021 

0.0023 

0.0018 

0.0014 

0.0014 

0.0014 

0.0014 

0.0014 

Table  2:  Fitting  the  Between-Individual  Model 


/ 

It 

standard  error (ji) 

1 

-0.9022 

0.084186 

2 

-0.8336 

0.068770 

3 

-0.7184 

0.052690 

4 

-0.6567 

0.047560 

Table  3:  95%  Confidence  Intervals  for  the  Difference  between  Expected  Log 
Deterioration  Rates  of  Two  Groups _ 


group  l  (vendor  and  depot) 

-  group  l1  (vendor  and  depot) 

(lower  limit,  upper  limit) 

1  (A  and  1) 

- 

2  (B  and  1) 

(-0.2990,  0.1618) 

1  (A  and  1) 

- 

3  (A  and  2) 

(-0.3943,  0.0267) 

1  (A  and  1) 

- 

4  (B  and  2) 

(-0.4504,  -0.0406) 

2  (B  and  1) 

- 

3  (A  and  2) 

(-0.2989,  0.0685) 

2  (A  and  1) 

- 

4  (B  and  2) 

(-0.3542,  0.0004) 

3  (A  and  2) 

- 

4  (B  and  2) 

(-0.2122,  0.0888) 

Table  4:  Predicted  Proportion  Defective  at  year  5  and  Estimated  Time  for 

p  —  0.5 _ 

P)  (%)  t,  (year) 

5.34  12.1 

6.12  11.3 

7.83  10.1 


9.02 


9.5 
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